High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

New traveling wave solutions of Drinefel’d–Sokolov–Wilson Equation using Tanh and Extended Tanh methods

Traveling wave solutions are obtained by using a relatively new technique which is called Tanh and extended Tanh method for Drinefel’d–Sokolov–Wilson Equations. Solution procedure and obtained results re-confirm the efficiency of the proposed scheme.     

High-resolution semi-discrete Hermite central-upwind scheme for multidimensional Hamilton–Jacobi equations

We introduce a high resolution fifth-order semi-discrete Hermite central-upwind scheme for multidimensional Hamilton–Jacobi equations. The numerical fluxes of the scheme are constructed by Hermite polynomials which can be obtained by using the short-time assignment of the first derivatives. The extensions of the proposed semi-discrete Hermite central-upwind scheme to multidimensional cases are straightforward. The accuracy, efficiency and stability properties of our schemes are finally demonstrated via a variety of numerical examples.     

Modified fully discretized projection method for the incompressible Navier–Stokes equations

The stability and convergence of a second-order fully discretized projection method for the incompressible Navier–Stokes equations is studied. In order to update the pressure field faster, modified fully discretized projection methods are proposed. It results in a nearly second-order method. This method sacrifices a little of accuracy, but it requires much less computations at each time step. It is very appropriate for actual computations. The comparison with other methods for the driven-cavity problem is presented.     

An improved algorithm for the shallow water equations model reduction: Dynamic Mode Decomposition vs POD

We propose an improved framework for dynamic mode decomposition (DMD) of 2‐D flows for problems originating from meteorology when a large time step acts like a filter in obtaining the significant Koopman modes, therefore, the classic DMD method is not effective. This study is motivated by the need to further clarify the connection between Koopman modes and proper orthogonal decomposition (POD) dynamic modes. We apply DMD and POD to derive reduced order models (ROM) of the shallow water equations. Key innovations for the DMD‐based ROM introduced in this paper are the use of the Moore–Penrose pseudoinverse in the DMD computation that produced an accurate result and a novel selection method for the DMD modes and associated amplitudes and Ritz values. A quantitative comparison of the spatial modes computed from the two decompositions is performed, and a rigorous error analysis for the ROM models obtained is presented. Copyright © 2015 John Wiley & Sons, Ltd.     

Evolving a three-world framework for solving algebraic equations in the light of what a student has met before

In this paper we consider data from a study in which students shift from linear to quadratic equations in ways that do not conform to established theoretical frameworks. In solving linear equations, the students did not exhibit the ‘didactic cut’ of Filloy and Rojano (1989) or the subtleties arising from conceiving an equation as a balance (Vlassis, 2002). Instead they used ‘procedural embodiments’, shifting terms around with added ‘rules’ to obtain the correct answer (Lima & Tall, 2008). Faced with quadratic equations, the students learn to apply the formula with little success. The interpretation of this data requires earlier theories to be seen within a more comprehensive framework that places them in an evolving context. We use the developing framework of three worlds of mathematics (Tall, 2004, Tall, 2013), based fundamentally on human perceptions and actions and their consequences, at each stage taking into account the experiences that students have ‘met-before’ (Lima and Tall, 2008, McGowen and Tall, 2010). These experiences may be supportive in new contexts, encouraging pleasurable generalization, or problematic, causing confusion and even mathematical anxiety. We consider how this framework explains and predicts the observed data, how it evolves from earlier theories, and how it gives insights that have both theoretical and practical consequences.     

Generalization of tangential trace spaces of H(curl,Ω) for curvilinear Lipschitz polyhedral domains Ω

For curvilinear Lipschitz polyhedral domains Ω, explicit characterizations of the tangential trace spaces of H ¹(Ω) are presented. These extend the original characterizations given by Buffa and Ciarlet that hold on Lipschitz polyhedral domains with plane faces. The tangential trace spaces of H ¹(Ω) are fundamental for the definition, analysis and intuitive understanding of the trace spaces of H ( curl ,Ω) and therefore, more general characterizations of the latter are obtained at the same time. Copyright © 2013 John Wiley & Sons, Ltd.     

On analyticity of the ‐Stokes semigroup for some non‐Helmholtz domains

Consider the Stokes equations in a sector‐like C ³ domain . It is shown that the Stokes operator generates an analytic semigroup in for . This includes domains where the ‐Helmholtz decomposition fails to hold. To show our result we interpolate results of the Stokes semigroup in and L ² by constructing a suitable non‐Helmholtz projection to solenoidal spaces.     

Error Estimates for a Class of Elliptic Optimal Control Problems

In this article, functional type a posteriori error estimates are presented for a certain class of optimal control problems with elliptic partial differential equation constraints. It is assumed that in the cost functional the state is measured in terms of the energy norm generated by the state equation. The functional a posteriori error estimates developed by Repin in the late 1990s are applied to estimate the cost function value from both sides without requiring the exact solution of the state equation. Moreover, a lower bound for the minimal cost functional value is derived. A meaningful error quantity coinciding with the gap between the cost functional values of an arbitrary admissible control and the optimal control is introduced. This error quantity can be estimated from both sides using the estimates for the cost functional value. The theoretical results are confirmed by numerical tests.